Soft Dynamics simulation: 
normal approach of deformable particles in a viscous fluid 
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Discrete simulation methods are efficient tools to investigate the complex behaviors of complex fluids made 
of either dry granular materials or dilute suspensions. By contrast, materials made of soft and/or concentrated 
units (emulsions, foams, vesicles, dense suspensions) can exhibit both significant elastic particle deflections 
(Hertz-like response) and strong viscous forces (squeezed liquid). We point out that the gap between two 
particles is then not determined solely by the positions of their centers, but rather exhibits its own dynamics. 
We provide a new discrete numerical method, named Soft Dynamics, to simulate the combined dynamics of 
particles and contacts. As an illustration, we present the results for the approach of two particles. We recover 
the scaling behaviors expected in three limits: the Stokes limit for very large gaps, the Poiseuille-lubricated 
limit for small gaps and even smaller surface deflections, and the Hertz limit for significant surface deflections. 
Larger scale simulations with this new method are a promising tool for investigating the collective behaviors of 
many complex materials. 
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PACS numbers: 02.70.Ns, 82.70.-y, 83.80.Iz 



Foams, emulsions and granular materials are made of inter- 
acting particles, respectively bubbles, droplets and grains, in a 
surrounding fluid. A great deal of research revealed their elas- 
tic, plastic and viscous characters QUI]. How are these behav- 
iors related to the particle interactions? This issue is central 
in the current debate about constitutive laws, dilatancy JH01, 
and shear banding JH 0, 0, HI HI . Discrete simulation methods 
are efficient tools for such an investigation, since the motion of 
each particle is derived from the interactions directly. Molec- 
ular Dynamics (MD) for dry granular fioll . Stokesian Dynam- 
ics (SD) for viscous suspensions ifTTIl and Bubble Model (BM) 
for foams IU2I1 implement various adjustable particle interac- 
tions and provide the corresponding macroscopic behaviors, 
with good agreement with experiments lfl3i [14 [l5ll . Hence, 
implementing relevant interaction models is crucial. Accurate 
models are already available for either deformable elastic-like 
grains without surrounding fluid (MD for dry granular), or 
non-deformable particles in a viscous fluid (SD for suspen- 
sion). By contrast, the combination of both particle defor- 
mation and viscous flow has not been fully described yet, al- 
though it is central in such materials as foams and emulsions. 
BM provided a first step in that direction. Like with MD for 
grains, BM particles can slightly overlap with an associated 
normal elastic-like repulsion. Besides, a tangential viscous 
force arises from the shearing of the contact film. However, 
BM does not include any viscous effects for normal displace- 
ment, and is therefore not fully suitable for the normal ap- 
proach or separation of particles, which are essential in Tl 
and other topological processes. 

In this paper, we address the description of normal inter- 
actions between two elastic-like particles in a viscous fluid, 
as depicted in Fig. Q] We show that such a system exhibits 
two distinct dynamics for the center-to-center and surface-to- 



surface distances. As this feature was ignored so far in exist- 
ing simulations, we here introduce a new discrete numerical 
method, named Soft Dynamics, and implement the simple case 
of two similar particles for illustration purposes. We recover 
the scaling behaviors expected in Stokes, Poiseuille-lubricated 
and Hertz limits. 
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FIG. 1: Sketch of normal particle interaction. A force between elastic 
surfaces (modulus E) is transmitted partly through the fluid (visosity 
r)) and partly through possible remote interaction (labeled W). It can 
deform them (deflection 8), which affects the flow. The consequence 
is that the center-to-center distance X and the gap h between both 
surfaces exhibit two distinct dynamics. 

Model of interaction. Let us describe a minimal model of 
interaction between two identical, elastic spheres with Young 
modulus E and radius R within a Newtonian fluid of viscosity 
tj . We neglect both particle and fluid inertia, thereby assuming 
the system to be in a Stokes regime. Furthermore, we focus 
on normal motion. An external force F ext induces a relative 
motion of particles, so that both their inter-centre distance X 
and the gap h evolve in time. Particle surfaces are then sub- 
mitted to the viscous force F vls due to the motion of the fluid 
in the gap, and also to a possible remote interaction F adh that 
is not transmitted through the fluid. Submitted to these two 
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forces, particle deform elastically (surface normal deflection 
8, counted as positive while in compression). The force bal- 
ance (with no inertia) and the geometry yield the following 
equations of the system (Fig.Q~|i: 



pext pela pvis _j_ padh 

X = 2R + h-28. 



(1) 

(2) 

According to Hertz iflfjll . the elastic deformation of a sphere 
due to a force acting on a region of radius a is mainly located 
within a volume of size a, and scales like S/a, Hence, the 
total elastic energy stored in the contact scales like EaS 2 and 
the associated force reads: 



F da (a,S) = aEaS. 



(3) 



where the geometrical constant a is equal to 4/9 for elastic 
solids. The Young modulus E ranges typically from 100 GPa 
for metals to 100 MPa for polymers or 1 MPa forbiomaterials. 
By contrast, Eq. (01 is not so accurate for bubbles and droplets 
since their bulk is not elastic. Nevertheless, due to the surface 
tension, <7, their surface deflects in an elastic-like manner. The 
effective Young modulus is much smaller than for solids since 
it scales like a/R and thus ranges from 1 Pa to 10 4 Pa (a « 
10 _2 Nm _1 , S«l cm to 1 /im). Although other parameters 
can play a role in the elastic-like properties of bubbles and 
droplets, such as the number of contacts lfl7[[l8ll . the elastic- 
like force mainly depends on the deflection 8 and size a of the 
interacting region. 

The viscous force exerted by the fluid on the particle sur- 
faces, due to their normal motion h, is well established in 
three limits. First, the Stokes limit holds for very large gaps 
(/? ^> R): the viscous force 6nrjRh is transmitted through a 
region of size a w R. Secondly, for smaller gaps (/; < R), 
and negligible surface deflections (|5| <§; h), the interaction 
can be regarded as a Poiseuille lubrication flow between two 

spheres. Hence, the viscous force 67r ^ fi h is mainly transmit- 
ted through a region of size a where the surface-to-surface 
distance lies between h and 2h, yielding a \ / 2Rh (Fig. [2] 
a). Thirdly, the Hertz limit corresponds to elastic deflections 
larger than the gap (|5| h). In this case, the force is trans- 
mitted mainly throug h a region of size a determined by Hertz 
theory, a sa y/2R\8\, and the contact can be regarded as two 

parallel disks, so that the viscous force reads ^-h. In or- 
der to interpolate between these three limits, we express the 
viscous force as: 



F vis (a,h,h) = -Ch, 



(4) 
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is the effective normal friction coeffi- 

2/7 



where £(a,h) 

cient of the contact, with h = Rh/(R + h). Other choices are 
possible for this coefficient, but it will at least depend on the 
contact geometry through h and a. We choose an expression 
for the size of the interacting region that interpolates between 
the asymptotic expressions valid in all three regimes (Stokes, 
Poiseuille and Hertz): 



a{h,8) = \2R(h+\8\) 



(5) 



Other choices for a are possible, but a will depend at least on 
h and 8. 
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FIG. 2: Scheme of the region through which both particles inter- 
act (bold lines), (a) In the Poiseuille lubricated limit for sligthly 
deflected surfaces (S <C h), the size a of the interacting region de- 
pends mainly on h and marginally on 8. (b) By Contrast, in the 
Hertz limit for significant deflections (S h), it depends mainly on 
S and marginally on h. 

A remote interaction between both particle surfaces can 
also be present. It may include steric repulsion (with a closest 
approach distance related to the size of the surface asperities 
in solid grains), surfactant repulsion or, more generally, dis- 
joining pressure in foams or emulsions lfj~9n . It may also in- 
clude some adhesion. Various models can be developped for 
such remote forces, which will at least depend on the gap h 
and on the size a of the interacting region: F adh (a,h). In the 
present work, for simplicity, we omit such an interaction. 

Let us now consider two particles with elastic deflection 
and viscous interactions. This situation evokes a Mawxell 
fluid whose elastic and viscous coefficients would depend on 
the gap h and on the size a of the interacting region. We can 
derive the expected scaling behaviors of 8, a and h with re- 
spect to h in several asymptotic regimes, and express them 
simply (see Fig.[3]a,b) in terms of a typical Stokes velocity V 
and a renormalized applied force K defined by: 



V = F ex '/(r]R), K = F ex '/(ER 2 ). 



(6) 



We do not address large elastic deformation (8 ~ R) that cor- 
responds to K > 1. In the Stokes regime (h ^ R) the de- 
flection 8 can be neglected in a first approximation, and we 
find h ~ V and a ~ R. According to Eq. (O, we then ob- 
tain the deflection 8 ~ K. The crossover to the Poiseuille- 
lubricated regime occurs when the gap h becomes compara- 
ble to R, while remaining much larger than 8. Again, in the 
Poiseuille regime, 8 can be neglected at first, and we have 
a ~ VRh- Thus, from Eqs. and ©, h ~ Vh/R. Then, 
from Eqs. (O and ©, we get 8 ~ \J (R? /Ii)k. The crossover 
between Poiseuille and Hertz regimes occurs when the elas- 
tic deflection becomes comparable to the gap 8 ~ h, which 
corresponds to K ~ (/i/7?) 3 / 2 . In the Hertz regime, a ~ \/~R8. 
According to Eqs. (0 and ©, this implies a ~ Rk i I^. Finally, 
using Eqs. © and ©, we obtain h ~ V K^ 4 ^(h/R) 3 . 

Soft Dynamics simulation. Let us now show that the gap /; 
is not only dictated by the center-to-center position X, but ex- 
hibits its own dynamics. The system evolution is determined 
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FIG. 3: Scaling behaviors of two elastic-like sphere in a viscous fluid, (a) Scaling expected from theory in the Stokes, the Poiseuille-lubricated 
and the Hertz limits; the elastic deflection S/R, the contact size a/R and the gap rate of change h/V are function of the current gap h/R and 
renoralized force K. (b) Phase-diagram expected from theory for the crossover (large lines) between the three limits ; lines at constant K" come 
from Soft Dynamics simulations (see Fig.0. (c-f) 8(h), h(h) and their power law also from these simulations. Theoretical phase-diagram is 
replotted on the 3D plot basis. 



by Eqs. ((TJ and ||2), and by the expressions for the forces 
F ela (a,8), F adh (a,h), F vis (a,h,h) and a(h,S). The viscous 
force (Eq. @]i and force balance (Eq. [U directly yield the gap 
evolution: 

pext padh 

h= . (7) 

This is a purely viscous behavior with an effective friction 
£(a,h) and an applied force F ext — F adh which corresponds to 
the part of the total force that is transmitted through the fluid. 
Obtaining the dynamics of X is also straightforward, although 
it requires more steps. To this aim, it proves helpful to express 
time derivative of a and F ela in terms of their partial deriva- 
tives (denoted by a subscript, for example FfJ" = dF ela /da): 

a = a h h + a s 8, F el " = F? a a + F e s la 8. (8) 

Substituting d and 25 =h — X (from Eq. |2]i into the expres- 
sion of F ela yields F el " = -tfX + h(^ + Ff a a h ), where = 
(F t f a a s +Fl u ')/2. Then, using the force balance F ela = F ext 
and the expression for h (Eq.[7j leads to the dynamics equation 
foxX: 

pext pext padh 



with 9 = c €t,l( c € + Ff a a h ). The above expression evokes 
a Maxwell behavior, involving the time derivative of the ex- 
ternal force, with stiffness coefficient c <§ and friction coeffi- 
cient It thus differs from the purely viscous behavior of the 



gap by itself (Eq. O. Both dynamics of X and h are here ex- 
pressed in a generic and explicit manner in terms of the partial 
derivatives Ff'| and a/, g and the effective friction coefficient 

£, and also a possible remote force F" dh . The expressions for 
F ela (a,8),F adk (a,h) anda(M) 

can therefore be changed at 
will to describe various physical contact behaviours. The only 
restriction is that the viscous force must be linear in h, consis- 
tently with the no-inertia assumption. 

Knowing the distinct dynamics of h and X, we can now pro- 
pose a method to simulate their time evolution in the presence 
of a known external force F ext (t) (see Fig.|4|i. Initial values of 
X and h yield the deflection 8 according to Eq. (0, and then 
the size of the contact a(h,8). Quantities £, ^ and Q> are ex- 
plicitely given by the expression of the forces and a. Finally, 
the external force F ext and its time derivative are needed to 
obtain both h and X. By integration, new values for h and X 
are obtained. 

We implemented the Soft Dynamics method in the simple 
case of two identical particles. They are initially separated 
by a large distance (h/R = 100), we then subject them to a 
constant, compressive force F ext , and the gap decreases con- 
tinually, h may even tend to zero asymptotically since in this 
case, for simplicity, we have not included the remote force 
padh p or var j ous 

applied forces (10~ 6 < K < 0.3), Fig. [3J) 
shows the range of gaps h reached throughout the approach. 
Plotting 8 and h as a function of the current value of the gap 
h allows us to extract the corresponding scaling behaviors. 
Figs. [3j;,d represent function 8(h) and the associated expo- 
nent X\(h) such that one locally has 8(h) h i. In both the 
Hertz and the Stokes regimes, the deflection 8 depends only 
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FIG. 4: Method for the Soft Dynamics simulation. From current 
center-to-center distance X and gap h, their respective time evolu- 
tions are explicitely calculated using the value of the external force 
F ext and its time derivative. The intermediate calculations involve 
analytical expressions that result from the physical description of the 
system (solid rectangular boxes). The dashed rectangular box (with 
and &) represents intermediate calculation steps. 



on the contact force, with the usual 2 /3 -power law in the Hertz 
regime and a linear law in the Stokes regime where the inter- 
acting region spans over one entire side of the particle. In the 
Poiseuille-lubricated regime, the width of the interacting re- 
gion depends on the gap h: it causes 8 to vary like {h/R)~ l l 2 . 
Function h(h) is represented on Figs. [3^,f together with the 
associated local exponent: h(h) h^ 2 . In the Stokes regime, 
the large gap h^$> R implies that the interaction, mediated by 
the surrounding liquid, does not depend on distance: h is then 
independent of h. By contrast, in both the Poiseuille and the 
Hertz regimes, the gap thickness affects the liquid flow, and /; 
varies respectively like hoc h l and h <x /i 3 . Note that the range 
of gaps h where k\ w — 1/2 and A2 ~ 1 shrinks and vanishes 
as the renormalized force k is increased: strongly compressed 
particles cross over from Stokes to Hertz regimes directly as 
h deacreases. All these observations are consistent with the 
theoretical predictions summarized on Fig.[3k. 



The conclusion of this work is that particles interacting in a 
viscous fluid can exhibit two distinct dynamics, which govern 
their overall motion on one hand and the gap evolution on the 
other hand. This behavior results from the interplay between 
the underlying Poiseuille and Hertz mechanisms: the viscous 
flow in the gap affects the transmitted force and hence the sur- 
face deflection, which in turn affects the flow in the gap. In 
Stokes and Poiseuille regimes, grain surfaces are remote and 
the deflection can be neglected: the gap dynamics is mainly 
dictated by the motions of the particle centers. However, as 
soon as the deflection becomes larger than the gap, both dy- 
namics must be distinguished. This is easily encountered with 
soft particles, but can also occur with rather hard particles as 
soon as the distance h becomes small enough. It seems that 
this feature had not been implemented so far: we here pro- 
vided a new simulation method, called Soft Dynamics, which 
expresses each of these two dynamics generically in terms of 
the forces involved (elastic, viscous as well as possible remote 
interaction). Here, we merely check that it describes the nor- 
mal approach of two particles in a viscous fluid correctly. But 
Soft Dynamics^ can be implemented for many multicontact- 
ing particles I20I1 . Hence, like MD for dry grains and SD for 
diluted suspensions, Soft Dynamics should be able to provide 
new insights into the behavior of materials made of soft, close- 
packed units. Bubbly liquids, vesicules, foams, emulsions, 
micellar solutions and dense suspensions could thus be inves- 
tigated, by including the specificities of their interactions. 
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